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Abstract 

We study two-component single-file diffusion inside a narrow channel that at its ends is open and 
connected with particle reservoirs. Using a two-species version of the symmetric simple exclusion 
process as a model, we propose a hydrodynamic description of the coarse-grained dynamics with 
a self-diffusion coefficient that is inversely proportional to the length of the channel. The theory 
predicts an unexpected nonequilibrium phase transition for the bulk particle density as the external 
total density gradient between the reservoirs is varied. The individual particle currents do not in 
general satisfy Fick's first law. These results are confirmed by extensive dynamical Monte-Carlo 
simulations for equal diffusivities of the two components. 
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I. INTRODUCTION 



One-dimensional exclusion processes belong to the most studied models in non- 
equilibrium statistical mechanics Their applications are manifold. Among others, 
the symmetric exclusion process (SEP) plays a role in diffusion where particles, confined in 
a narrow tube, are not allowed to pass each other Q|. This kind of diffusive restriction is 
referred to as single-file diffusion and differs qualitatively from normal diffusion described 
by Fick's law. Whereas in the latter case the mean-square displacement of a single particle 
grows proportional to time, diffusion is much slower in the single-file case due to mutual 
blocking of the particles. The mean-square displacement grows (for late times) proportional 
to the square root of time. The anomalous behaviour of the mean-square displacement 
usually serves as an experimental indication for the occurrence of single-file diffusion. This 
requires to trace a single or more particles which implies to label a certain subset of particles 
without changing the diffusion properties. This corresponds to having a two-species particle 
system with identical diffusion coefficients [J]. Single-file diffusion is a generic phenomenon 
observed many years ago for molecules diffusing in the channels of certain zeolites More 
recently, single-file behaviour has been demonstrated in the transport of colloidal spheres 
confined in one-dimensional channels jf|. Moreover confined ID random motion plays a role 
in narrow carbon nano tubes, in biological systems like molecular motors or in non-physical 
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i as automobile traffic flow [zj|. Also the famous reptoii model by Rubinstein 



Hi] for the motion of single polymer chains, is a lattice gas model of this 
kind. Further motivation for employing single-file diffusion with multiple species comes from 
recent two-species measurements in zeolites jt^]. Here, a mixture of toluene and propane 
was adsorbed into different zeolites. The authors measured the temperature dependent out- 
flow and noticed a trapping effect, i.e. in a couple of zeolites the stronger adsorbed toluene 
molecules influence and control the outflow of propane. 

In [13I the authors review the Maxwell-Stefan theory describing the diffusive behaviour 
of a binary fluid mixture where the total current, i.e. the sum of both species, is zero. 
The particle-particle interaction is taken into account by including a friction between the 
species being proportional to the differences in the velocities. This approach does not apply 
to single-file diffusion in a finite system where, as shown below in the framework of the 
symmetric simple exclusion process (SEP), the self-diffusivity of single particles plays an 
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FIG. 1: Three-state symmetric exclusion model with open boundaries. 

important role in the description of the macroscopic behaviour. 

The SEP with one species of particles where classical particles with hard-core repulsion 
diffuse on a finite lattice is well understood At both ends the chain 



is connected to a particle reservoir. One is interested in stationary-state properties like 
the density profile determined by the reservoir or the stationary particle current as well as 
the time evolution of the particle density and relaxation towards the stationary state. Our 
approach for an adequate description of the two-species SEP is a master equation description 
from which we derive an ansatz of coupled partial differential equations for the macroscopic 
density profile. 

We consider a one-dimensional lattice with L lattice sites (Fig. [I]). Each site i can be 
empty or occupied by a particle of type A or B. Due to hard-core interaction any site 
carries at most one particle. Particles can hop to nearest neighbour sites (provided the 
target site is empty) according to the constant hopping rates Da/b- Hence, Da (Db) is the 
probability of an A (B) particle to attempt a jump per unit time. The model is defined by 
random sequential update which forbids simultaneous hopping events. Let (foj) count the 
A (B) particles on site i. Then the densities are the expectation values of the respective 
counters: (aj) = Pa^), (h) = Psij)- The probability of finding no particle at site i is 
(vi) = (1 — ai — bi). When we consider a chain with open boundary conditions, particles are 
injected and removed according to the boundary rates oa/b, 1a/b, Pa/b and 5a/b, following 
the notation of The attempt probability per unit time for an A particle to enter the 
system at the left boundary is a a- It leaves the channel at the left boundary according to 
7a as illustrated in Fig. [D The other boundary rates are defined similarly. 

By writing this process in terms of a quantum Hamiltonian formalism Py, the system 
evolves in time according to the master equation 



\P(t) >= -H\P{t) > (1) 
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with the generator 



L-l 



H = bi + b L + ^ hi,i+i- 



(2) 



i=i 



For an explicit representation of the generator we denote the state of a given site % by the 
three basis vectors 

|0 >= |0 >= 



\A >= II >= 








1 



V / 



B >E 



1 >= 








v 1 / 



(3) 



corresponding to having an A, no particle or B at site i, respectively. Let E%' y be the 3x3 
matrix with one element located at column x and row y equal to one. All other elements 
are zero: (E^ ,v ) a ^ = 5 Xja 5 y ^. The operator for annihilation (creation) of an A particle at 
site k is = El' 2 (a£ = E 2 ' 1 ) and for annihilation (creation) of B is = E^ 2 (b£ = E^ 3 ). 
Finally, the number operators are a k = E^ 1 , b k = E^ 3 , v k = 1 — a k — b k . This allows to 
compose the generator of the process. In this representation the boundary matrices are: 



&i = a A {vi - at) + a B (vi - bf) + 7a(oi - a ± ) + j B (bi - b ± ) 
h = 5a{v l - at) + 5 B (v L - &£) + j3 A {a L - a£) + j3 B {b L - &£) 

Hopping in the bulk between site i and i + 1 occurs according to 



(4) 
(5) 



h iji+ i = D A (aiV i+ i + Via i+ i - a i af +1 - afa i+1 ) 

+ D B (biV i+1 + Vib i+1 - bfbr +1 - b~b+ +1 ). (6) 

The model is now well defined. Let us proceed by deriving some equilibrium properties of 
the process. 



II. EQUILIBRIUM PROPERTIES 

The open system allows for particle exchange at the boundaries. The system is ergodic 
and will relax to a unique stationary state determined by the boundary rates. The stationary 
state \P* > does not evolve in time and must obey 

H\P* >= 0. (7) 
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FIG. 2: Three-state symmetric exclusion model with open boundaries - reservoir picture. 



Let us seek a product ansatz for the equilibrium state of the form 

(8) 



P* >= 



' a * 
1 



a + b) 1 



The normalization factor of jEJ) ensures conservation of probability, i.e. ensures that the 
probability of finding the system in any state is one. Plugging the ansatz ([Sj) into ((3) 
provides a set of equations for the boundary rates and one finds a= ^ = ]^ib=^ = j^. 
Taking into account the normalization determines the A and B particle equilibrium densities 

n _ Ja _ Pa /q\ 

PA I + + £s i , 6a , &b_ Ky > 

ia ^ 7s 1 ^ Pa ^ Pb 

Pb = , 2a T «b = i , ii. , fe- ( 10 ) 

7A ^ 7s 1 ^ /3a ^ Pb 

Besides giving the bulk equilibrium densities, the two equations above provide a recipe 
of how to translate the picture of inserting and deleting particles on boundary sites into 
a picture of constant reservoirs at the ends, Fig. [2l In equilibrium the bulk contains no 
correlations between different lattice sites and the same holds for the boundary and their 
adjacent sites. Therefore, jumping from a boundary site into the chain occurs proportional 
to the respective hopping rate and proportional to the single species boundary density. 
Therefore, given a set of constant boundary rates, Eqs. ([9]) and (fTOl) define the densities of 
a virtual particle reservoir at the respective boundaries. 

Note that for this interpretation the left reservoir densities p A , p^ do not need to be 
equal to their fellows on the right edge (p A , p^). In this case the system evolves towards a 
correlated non-equilibrium stationary state with non- vanishing particle currents. The second 
and last terms of ([9]) and ( TTOl) are then stationary densities on the left and right edge of the 
system. This parameterisation satisfies (El) and ( fTOl ) if 



a A/B = D A/B p A/B , 7 A/B = D A/B (l - p A - p B ) (11) 
S A/B = D A/B p + A/W f3 A/B = D A/B (1 - p\ - p%). (12) 
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III. HYDRODYNAMIC LIMIT 



The average densities (aA and (bi) satisfy the equations of motion 4 (cij) = — (aiH), 
Pi 

^ (bi) = — (hH) (cf. [2J). This provides the Master equations for a single site, 

(aA = L> A + (a i+1 Vi) - (tiiUi+i) - (aifi-i)) (13) 

= £> s ({fti-iUi) + - - (6iUi_i)) • (14) 

From now on we discuss the case of equal hopping rates Da = Db = D. Imagining the 
individual particle species A and B to be distinguishable by a "colour" (in an abstract sense) 
we shall refer to the total particle density (averaged over A and B particles) as colourblind 
density. 

Let us first assume an infinite system and do not care about boundary sites. But still, 
in this form the equations of motion are not integrable. Replacing the joint probabilities by 
products of expectation values, according to a mean field ansatz which has been proven to 
be useful in other systems, fails. However, an exact equation containing no correlators can 
be achieved from a sum of both 

j t ((a,) + (k)) = D ((a^> + (a i+1 ) - 2 (a t ) + (6<_i) + (b i+1 ) - 2 (b, t )) . (15) 

The right-hand side contains a second-order difference for both species individually. (TT5T) 
is the discrete analogue of the diffusion equation for the colourblind macroscopic profile. 
Introducing a lattice constant a and replacing i by the continuous variable x — -, transforms 
f fT5"l) for the hydrodynamic limit of vanishing lattice constant a into 

d t (p A (x, t) + p B (x, t)) = Ddl(p A (x, t) + p B (x, t)). (16) 

for the macroscopic particle densities pa(x, t), pB(x,t). 

Following the argument of [3| we make an ansatz for the dynamics of a single particle 
localized at position x. For a short-time region this particle acts as a tracer particle in the 
background of other particles with the self-diffusion coefficient Ds{x,t). Going beyond [3| 
we argue that for the finite-size problem with open boundaries Ds is given by expression 

D s = D l - pA - pB l -. (17) 

PA+ PB L 

derived originally for a periodic lattice [16J. Additionally, the test particle is subjected to a 
drift b caused by the evolution of the entire system towards its stationary state. For a good 
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intermixed background one would expect the drift velocity b to be the same for both species 
of particles. We thus arrive at the ansatz 

d t p A (x, t) = d 2 x D s p A (x, t) - d x bp A (x, t) (18) 
dtPB(x, t) = d 2 x D s p B (x, t) - d x bp B (x, t). (19) 

The self-diffusion coefficient D$ as well as the drift b are functions of p A and pB and hence, 
depend implicitly on x and t. The drift term can be determined by using the colourblind 
exact result (TT6l) and one finds 

b = ~d x \p(D s - D)} . (20) 

This completes the hydrodynamic description of the two-component symmetric exclusion 
process with open boundaries that we propose. A derivation of Dq on a finite lattice with 

n n 

two different particle species [20j will be presented in a forthcoming paper [21|. 

IV. STATIONARY STATE 

It is a significant property of particle systems with open boundaries that they can relax to 
a steady state with non- vanishing particle currents. The stationary state of the colour-blind 
profile p = p A + ps is linear with the slope being determined by the sum of the boundary 
densities on the left (p~ = P a + Pb) and on the right edge (p + = p\ + p^) of the system 
which is manifest by (fTBl) . 

P = P- + (P + " P~) j- (21) 
Integrating ([TBI once for vanishing time derivative yields 

where j A is the constant A-particle current. Absorbing p~ A into the integration constant 
yields the solution 



Pa{x) = — 



L 3AP ( - Lj A p \ fl - p{x) 

~Zs t I P A + 



P + ^p-. (23) 



D{p+-p-) \ D(p+-p-)J V l-p- 
The first term is linear in x and describes the bulk region. The nonlinear second term 



describes a boundary layer, first observed numerically in the Rubinstein-Duke model [U| 



for different boundary rates. Our analysis shows that the length of the boundary layer does 
not scale with system size. This can be seen by rewriting ( l23l) for sufficiently large L and 
assuming p + > p~: 



P + ~ P x \ ( l 3aP , ( - , Lj A p 



The localization length 

{ = (25) 

does not depend on the system size. In the limit of infinite L the relative size of the 
boundary layer vanishes and the linear solution connects to the reservoir densities by a 
jump discontinuity at one of the edges. For p + > p~ the exponential in (124|) dominates for 
small x and the discontinuity is located on left boundary. The case p + < p~ is similar, but 
the discontinuity is at the right edge. 

The case of equal reservoir densities p + = p_ = p has to be treated separately. The 
self-diffusion coefficient is constant, hence, the b in ([TBI vanishes. Integrating ( fT8l ) with 
Pa(0) = yields the linear density profile 

x 

Pa{x) = P A + (p\ - Pa) ^ (26) 

and a similar expression for the density of 5-particles. 

Fig. [3] shows the A and B particle densities obtained from Monte Carlo simulations 
(symbols) and the theoretical stationary state solutions (solid lines). The explicit expressions 

reservoir densities 



111 ]. This choice is 



for the particle currents are given below. We apply the same set o 

Pa = Pb = V3, Pb = i +e -i +e and Pa = i+l-i +e used in Fi S- 9 of 
motivated by boundary rates used in the Rubinstein-Duke model for describing the tensile 

force acting at the chain ends of reptating polymers. The different lattice sizes a) L = 50 

and b) L = 200 demonstrate the finite-size character of the boundary layer. The theoretical 

solution does not contain an inflection point and deviates slightly from simulations in the 

immediate vicinity of the boundary. Nevertheless, an interesting observation captured by 

the theoretical description is confirmed. The solution has a minimum in the density profile 

(although not very pronounced in the sample A-profile of Fig. [3]) and, hence, as in the 

Rubinstein-Duke model there exists a region where one of the particle currents does not 

follow the direction of the density gradient. 
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FIG. 3: Stationary state with finite slope of the colour-blind density, a) L = 50, b) 
L = 100. The boundary densities are Pa = Pb = 1/3, Pa = 1+e -i +e and p\ 



We conclude by analyzing the behaviour of the current and the mean particle density in 
the system. Using (1231) . (l26l) and taking into account the A particle reservoir density on the 
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right edge gives the current 



3a 



P A p A ( 1-P 

d(p+-p~) 7^ p- yi-p 

i-r ^ 



.i-p- 



(27) 



1) 



V>pO>a-Pa)0--p) P + = P 
This has an interesting consequence. Considering large L (l27j) simplifies asymptotically to 

p+ < p- (28) 



£> (p + ~P )p A 



Hence, provided a finite slope of the colour-blind profile, the individual particle currents 
are proportional to 1/L, as for the single-component case, whereas for p + = p~ the currents 
vanish proportional to 1/L 2 . 



We make an interesting observation if the relation p~ A p\ = PaPb * s satisfied. For this 
particular case the individual density profiles are linear and the particle currents are just 
proportional to the respective density gradients (Fick's law). If the relation does not apply 
we observe a boundary layer inside which the current flows against the local gradient. Here 
Fick's law is violated. 

Finally, the mean A-density in the channel can be obtained by integrating (1241) . Asymp- 
totically for large L one finds from (l28l) 

> P~ 

< p~ (29) 
= P~- 

The mean A-density evaluated as a function of the boundary densities may have a disconti- 
nuity. Assume p~ and p~ A ^ p\ be fixed. When taking the limit p + — > p~ coming from small 
p + the total density approaches p~A — > p A - Taking the limit from the other site, ~p~A — > p\- 
There is a jump of the mean A-density when the colour-blind boundary densities become 
equal. Since the colour-blind density is the sum of A and B densities, this implies a jump dis- 
continuity also in the mean 5-density. Therefore there is a first-order nonequilibrium phase 
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2p+ P 
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Pa+Pa 
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transition in this boundary-driven lattice gas model for two-component single-file diffusion. 
Such a transition is not known for boundary-driven one-component systems. 
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